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molecular  process  responsible.  The  associated  jump  time  can  be  Inferred  only 
on  the  basis  of  a  specific  model  for  the  motion.  By  a  reasonable  estimate  of 
the  relation  between  T.  and  the  jump  time  T,  we  characterize  the  motion 
responsible  for  the  high  temperature  relaxation  by  the  jump  time. 

T  ~  2  x  10~ 16  exp  (120  000/RT) 
where  T  is  in  seconds  and  R  in  units  of  J  mol”*  deg”*. 

The  motion  cannot  be  unambiguously  identified  with  those  results :  rotations 
or  180  flips  of  the  amino  protons  are  possible  as  well  as  rotation  about 
the  TATB  three-fold  axis.  By  comparison  with  other  molecules,  we  can,  however 
rule  out  translational  diffusion  as  the  motion  contributing  to  this  relaxation 
above  150°C. 

The  proton  NMR  spectrum  shows  the  characteristic  Pake  doublet  due  to  the 
nearly  isolated  proton  pairs  of  the  amino  groups.  The  observed  dipolar 
splitting  predicts  an  "average"  amino  proton  separation  of  0.185  +  0.001  nm, 
in  contrast  to  previous  x-ray  results  of  0.1553,  0.1577  and  0.1726  +  0.02  nm 
for  the  three  distinct  pairs  in  TATB. 
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PROTON  RELAXATION  IN  1, 3, 5-TRIAMINO-2, 4, 6-TRINITROBENZENE  (TATB) 


I .  INTRODUCTION 

The  compound  1,3,5  triamino-2, 4, 6- trinitrobenzene  (TATB)  is  simple 
yet  intractable;  it  is  virtually  insoluble  and  sublimes  with  decom¬ 
position  at  elevated  temperature.  It  is,  however,  extremely  attractive 
in  its  role  as  an  insensitive  high  explosive.  The  origin  of  this 
insensitivity  is  not  presently  known. 

Compressed  powders  of  pure  TATB  and  also  plastic  bonded  TATB 
materials  are  found  to  expand  irreversibly  when  the  temperature  is 
cycled;  it  is  suggested  (1)  that  this  hysteresis  stems  from  stresses 
induced  by  the  anisotropic  volume  expansion  of  TATB.  The  mechanism 
for  the  large  c-axis  thermal  expansion  coefficient,  over  eight  times  that 
of  graphite,  is  not  understood. 

As  a  step  towards  identifying  the  character  of  the  molecular 
motions  which  may  relate  to  mechanical  insensitivity  and  thermal 
expansion,  we  have  performed  a  series  of  proton  relaxation  studies 
in  (solid)  TATB  over  the  temperature  range  of  approximately  -80  to 
240  C.  In  brief,  we  find  that  there  is  very  little  (detectable) 
molecular  motion  below  about  150°C.  From  about  150  to  200°C  we  observe 
a  motion  with  an  activation  enthalpy  of  120  kJ/mol  and  associated 
frequency^ prefactor  of  0.5xl0x  Hz.  At  205  C,  this  corresponds  to  a 
jump  time  of  around  5  ms,  while  at  room  temperature  five  days 
elapses  between  those  jumps.  We  find  this  activation  enthalpy  is  too 
small  for  translational  diffusion;  the  motion  may  involve  large  angle 
libration  or  rotation  of  the  amino  group  or  in-plane  rotation  of  the 
entire  molecule,  though  other  mechanisms  are  not  ruled  out  by  these 
data.  In  particular,  we  observe  no  evidence  for  nonsubtle  changes 
either  in  the  crystal  structure  or  in  the  nature  of  the  molecular 
motions  in  the  temperature  range  around  30°C.  Relaxation  data  at 
elevated  temperatures  up  to  the  decomposition  point  (350°C)  might  help 
to  distinguish  further  among  the  possible  models  for  motion  in  TATB. 

II.  EXPERIMENTAL 

TATB  powder  (micronized,  chlorine-free)  was  obtained  from  J.  R. 

Kolb  of  the  Lawrence  Livermore  Laboratory.  In  the  as-received  material, 
a  narrow  proton  NMR  line  near  the  center  of  the  Pake  doublet  is 
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ascribed  to  liquid  water  (see  Section  III).  This  water  was  removed 
by  heating  the  samples  to  200°C  for  0.5  to  2  hours  under  vacuum;  a  0.1% 
weight  loss  occurred. 

For  the  T2,  and  T^relaxation  data,  the  pumped  powder  was  then 
compressed  into  pellets  (0.55  cm  dia  x  1.2  cm)  under  about  14MPa  (2000  psi) 
pressure.  Samples  fit  into  glass  tubes;  a  copper-cons tantan  thermocouple 
junction  was  suspended  about  1  cm  above  the  pellet  and  the  sample  tube  was 
nominally  sealed  off  by  a  stopper. 

For  technical  reasons,  a  separate  probe  was  used  for  the  few  T, 

Ip 

measurements  reported  here.  In  this  case  a  separate  pumped  specimen  was 
pressed  into  a  1.25  cm  dia  x  3  cm  cylinder.  Thermocouples  were  placed  in 
access  ports  of  the  heavy  copper  probe  body;  subsequent  bench  tests  showed 
these  reflected  the  temperature  within  the  sample  coil  to  better  than  3°C. 
Few  T  data  are  presented;  the  sample  degraded  at  about  300°C,  and  damaged 
the  probe. 

A  Bruker  SXP  NMR  spectrometer  was  used.  For  the  T2,  ^  and  T^D 

measurements,  the  proton  Larmor  frequency  was  60  MHz  and  the  rf  field 
strength  was  150-200  kHz.  (Field  strengths  are  reported  in  frequency 
units  by  the  relation  f=*yHB/27r,  where  B  is  the  field  and  y^  the  proton 

magnetogyric  ratio.)  A  special  probe  was  required  for  the  T.  measure¬ 
ments;  this  operated  at  25  MHz  and  a  spin- lock  rf  pulse  amplitude  of 
50  kHz  was  selected. 

The  T2  data  were  determined  by  the  zero  crossing  of  the  free 

induction  decay  (fid) .  To  minimize  receiver  recovery  time,  the  fid 
was  observed  immediately  following  a  spin-locking  sequence  consisting 
of  a  7r/2  pulse  and  a  lOOys  phase-shifted  spin-lock  pulse.  The  phase  of 
the  initial  tt/2  pulse  was  alternately  changed  by  180°  while  the  signal 
was  synchronously  added  and  subtracted  to  computer  memory;  this  was  done 
to  discriminate  against  instrumental  transients  and  in  favor  of  the 
nuclear  signal.  The  start  of  the  fid  was  considered  to  occur  at  the 
very  end  of  the  spin- lock  pulse. 

The  solid  echo  pulse  sequence,  ir/2(0°)-T-ir/2(90o)-T-echo,  was  not 
satisfactory  for  T2  determination.  (In  this  notation,  the  figures  in 

parenthesis  give  the  relative  rf  phase.)  Though  the  echo  refocussing 
at  the  time  2t  is  perfect  for  an  isolated  spin  pair  (2,3),  the  amino 
proton  pairs  in  TATB  are  not  completely  isolated.  Many  of  the  con¬ 
tributions  to  the  fid  time  evolution  which  arise  from  outside  the  spin 
pair  will  lead  to  a  systematic  error  in  the  fid  zero  crossing  (3) . 

The  dipolar  spin-lattice  relaxation  time  T-D  was  determined  bg 
the  two  guise  method  of  Jeener  and  Broekaert  (4;,  'n,/2(0°)-T-TT/4(90  )- 
t-rr/4(90D).  The  separation  T  between  the  first  two  preparation  pulses 
is  to  coincide  with  the  steepest  part  of  the  fid;  T  was  set  to  12  ys, 
as  measured  between  the  rising  edges  of  the  pulses.  As  above,  the  phase 
of  the  first  pulse  was  alternated  from  0°  to  180°  to  reduce  transients. 


« 


In  TATB,  the  proton  Zeeman  spin-lattice  (longitudinal)  relaxation 
time  ^  is  quite  long,  of  the  order  of  30  sec.  In  the  T2,  and 

experiments,  the  magnetization  was  not  allowed  to  return  to  its  thermal 
equilibrium  value.  Instead,  to  avoid  any  ambiguity  about  what  magneti¬ 
zation  was  inspected  by  the  sequence,  each  sequence  began  by  destroying 
all  magnetization  along  the  Bq  (z)  axis  by  a  comb  of  12  ir/2  pulses 
spaced  80  ysec  apart.  Then  a°10  sec  delay  was  allowed  for  the  magneti¬ 
zation  to  recover. 

The  T^  determination  proceeded  also  by  destroying  magnetization  by 

the  12  pulse  comb  and  then  monitoring  the  magnetization  recovery  of  a 
two-pulse  solid  echo. 

For  T^p»  the  decay  of  the  fid  was  monitored  for  varying  spin-lock 

pulse  lengths,  at  constant  rf  field  of  50  kHz.  The  rf  field  depend¬ 
ence  was  not  studied  here. 

III.  RESULTS  AND  DISCUSSION 

A.  Spin-spin  Relaxation  (T2) 


The  proton  NMR  spectrum  of  TATB  is  shown  in  Fig.  1;  the  narrow 
water  peak  in  the  as-received  material  (a)  has  essentially  disappeared 
after  heating  at  200  C  under  vacuum  (b)  with  a  concomitant  weight  loss 
of  0.1%.  In  Fig.  1  the  linewidth  of  the  narrow  resonance  Is  artifically 
broad,  due  to  truncation  of  the  free  induction  decay  (fid).  Other 
spectra  (not  shown)  obtained  with  a  longer  acquisition  time  give  a  more 
realistic  value  for  the  linewidth  of  700+50  Hz.  This  is  substantially 
greater  than  magnet  inhomogenity  (40Hz)  or  bulk  susceptibility  effects. 

We  infer  that  the  narrow  resonance  arises  from  water  which  is  partially 
constrained  by  the  TATB  powder.  The  water  peak  at  6=4. 1+0. 8  ppm 
(relative  to  TMS)  is  essentially  unshifted  from  neat  water  at  6=4.7  ppm. 

The  spectrum  of  Fig.  lb  represents  the  classical  and  character¬ 
istic  Pake  doublet  (5) ,  here  arising  from  the  dipolar  coupling  of  the  amino 
proton  pair;  the  spectrum  is  broadened  slightly  by  mechanisms  discussed 
below.  The  center  of  the  spectrum  is  shifted  approximately  34+8  ppm 
downfield  from  TMS.  This  shift,  also  observed  in  TATB  by  Fukushima  (6), 
reflects  the  substantial  hydrogen  bonding  in  the  crystal. 

We  look  for  any  change  with  temperature  in  the  dipolar  coupling; 
such  a  change  indicates  the  presence  of  a  molecular  motion  modulating 
the  coupling  at  a  frequency  greater  than  or  equal  to  the  coupling  frequency 
(about  30kHz  for  the  dominant  splitting  in  the  Pake  doublet).  Now,  Fig.  1 
is  not  a  perfectly  faithful  rendition  of  the  TATB  spectrum;  all  such  spectra 
are  distorted  by  imperfections  arising  from  the  overloading  of  the  NMR 
receiver  following  the  end  of  the  rf  transmitter  pulse.  Here  that  overload 
time  is  6.5  ys.  As  an  alternative  to  comparing  details  of  distorted  spectra, 
it  is  more  convenient  instead  to  characterize  the  entire  spectrum  by  the  first 
zero  crossing  of  the  fid  and  also  by  the  average  spectral  broadening 
produced  by  dipolar  coupling  outside  the  amino  group.  For  a  single 
Pake  doublet  (e.g.  for  a  powder  of  proton  pairs,  each  with  identical 
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intra-pair  distance) ,  the  time  dependence  of  the  normalized  free 
induction  decay  is  (7) 


G(t)  *  exp  [  -  1/2  g2t2]  f (at) , 


U] 


where 


f (at)  =  a(8Tr)1/2  ( cosat  C  (z1/2)  +  sinat  S(z1/2)J.  [2] 


1/2, 


1/2, 


(at) 


1/2 


(at) 


1/2 


Here  (4/5)  a  is  the  second  moment  of  the  intra-pair  coupling  and 

a  =  3y2  fl/(4r3) ,  [3] 


where  r  is  the  pair  separation.  The  second  moment  contribution  of 
all  the  other  spins  (the  extra-pair  contribution)  is  8  and  bjr  assumption 
its  contribution  to  the  fid  is  Gaussian  as  in  Eq.  [1].  The  argument  of 
the  Fresnel  sine  and  cosine  integrals  S  and  C  is  z  =6at/iT. 

The  advantage  of  the  zero  crossing  measurement  (G(t  )=0)  is 
easily  seen  from  Eq.  [1]:  the  zeros  of  G(t)  are  determined  completely 
by  the  intra-pair  coupling  a.  Further,  the  first  zero  crossing  occurs 
at  tQ*2. 160/a;  the  experimental  value  for  TATB  is  tQ^  25ys  and  the 

zero  crossing  event  is  well-removed  from  the  receiver  recovery  period; 
thus  tQ  is  less  subject  to  error  than  the  splitting  in  the  Pake  spectrum. 
Once  a  is  experimentally  determined,  the  extra-pair  contributions  (8  ) 
can  then  be  determined  by  dividing  G(t)  by  f(at). 

For  a  pumped,  compressed  powder  of  TATB  we  have  determined  the 
first  zero  crossing  with  its  associated  intra-  pair  second  moment  as 
a  function  of  temperature.  The  second  moment  results  are  shown  in 
Fig.  2  and  Table  I  while  zero  crossing  data  is  also  collected  in  Fig. 

5.  The  intra-pair  second  moment  ranges  from  8  to  10  G  while  the 
extra-  pair  contribution  is  about  3G  .  The  second  moments  are 
essentially  temperature  independent  over  the  range  -80  to  230  C, 
indicating  that  in  this  temperature  range  there  is  no  motion  which  modulates 
the  dipolar  couplings  at  a  rate  greater  or  equal  to  the  dipolar  coupling 
frequency  of  about  30  kHz.  (For  example,  if  the  TATB  molecule  were 
rotating  rapidly  about  its  three-fold  axis,  the  intra-pair  second 
moment  would  drop  by  a  factor  of  four,  as  would  any  other  intramolecular 
contributions) .  While  there  is  a  substantial  thermal  expansion  along  the 
c-axis  in  TATB,  the  proton  pairs  are  essentially  in  the  ab  plane  and 
their  separation  is  less  sensitive  to  thermal  expansion.  Using  the 
position  vectors  of  Cady  and  Larson  (8)  and  the  temperature  dependence 
of  the  unit  cell  parameters  of  Kolb  and  Risso  (1),  we  have  calculated 
that  there  is  less  than  a  2%  change  expected  in  the  intra-pair  second 
moment  over  the  temperature  range  -26  to  230°C.  Although  no  x-ray  data 
are  available  for  104  to  230  C,  in  this  range  we  used  the  average  expansion 
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coefficients  over  -54  to  104°C  together  with  the  values  of  the  unit 
cell  angles  for  104  C. 

From  the  experimental  value  of  the  intra-pair  second  moment,  we 
can  infer  the  pair  separation  via  Eq.  [3].  At  26°C  the  zero  crossing 
time  is  2422+0. 5]is  corresponding  to  an  intra-pair  second  moment  of 
8.9  +  0.4G  ;  the  experimental  value  for  the  extra-pair  moment  is 
3. 1+0. 5  G  ,  calculated  by  the  route  just  described.  The  pair  separ¬ 
ation  is  then  0.185  +  0.001  nrn.  This  value,  which  presumes  a  single  value 
for  all  the  pair  separations,  may  be  compared  to  the  x-ray  results  of 
Cady  and  Larson  (8) .  There  the  proton  separations  in  the  three  amino 
groups  are  0.1553,  0.1577  and  0.1726  +  0.02  nm,  where  we  have  calculated  the 
separations  from  the  atomic  position  vectors  given  by  the  spherical 
refinement  of  Ref.  8.  The  proton  distances  are  rather  poorly  identified 
by  the  x-ray  data  and  from  the  quoted  errors  in  the  position  vectors 
(8)  we  infer  that  the  associated  error  in  the  pair  separation  is 
approximately  +0.02  nm,  for  each  pair,  as  indicated.  Intra-pair  second 
moments  based  on  the  x-ray  determined  separations  are  shown  in  Table  I. 

The  NMR  estimate  of  0.185  +  0.001  nm  is  consistent  with  one  of  the  x-ray 
values  and  slightly  outside  the  range  of  uncertainty  for  the  two  other 
pairs. 


As  the  uncertainty  in  the  x-ray  determined  proton  separations  is 
so  large,  it  is  unclear  that  there  are  actually  three  widely  different 
pair  separations  and  so  no  attempt  was  made  here  to  scrutinize  the 
NMR  spectra  for  the  presence  of  three  distinct  pairs.  If  the  pair 
separation  were  not  unique,  then  the  zero  crossing  time  would  represent 
an  intermediate  average  over  the  three  pair  contributions  to  the  fid. 

As  all  of  the  x-ray  distances  fall  below  the  NMR  determination,  we  infer 
that  some  if  not  all  the  x-ray  proton  pair  separations  are  too  short. 


We  have  calculated  the  extra-pair  proton-proton  second  moment  from 
the  x-ray  data:  by  definition  this  is  the  second  moment  contribution  to 
each  of  the  six  distinct  TATB  amino  protons  made  by  all  other  protons 
in  the  crystal  but  excludes  the  proton  within  each  respective  pair. 

The  average  over  the  six  protons  is  just  the  extra-pair  contribution  to 
the  total  proton  second  moment.  The  contributions  of  near  neighbors 
out  to  1.01  nm  were  considered  and  we  estimate  this  truncation  under¬ 
represents  the  actual  second  moment  by  no  more  than  3%.  The  extra-pair 
second  moments  are  presented  in  Table  I  for  protons  numbered  1-6  in  the 
notation  of  Ref.  8.  These  moments  are  more  reliable  than  the  intra-pair 
moments,  as  they  involve  the  well-determined  unit  cell  dimensions  in 
addition  to  less-determined  intramolecular  proton  distances.  The  average 
extra-pair  interaction  is  1.50  G  which  is  substantially  less  than  the 
2 

3G  measured  experimentally  (Fig.  2). 


However,  the  dipolar  field  of  the  amino  nitrog^g  m^y  also  contribute 
to  the  proton  second  moment.  The  magnitude  of  the  H  coupling  is 

difficult  to  assign  for  a  number  of  reasons.  As  the  N  nucleus  has 
a  quadrupole  moment  and  sees  a  very  unsymmetrical  electronic  environment 
in  TATB,  the  quadrupolar  splitting  improbably  of  the  order  of  l-3MHz, 
which  is  not  too  different  than  the  N  Larmor  frequency  of  4.3  MHz  at 
this  magnetic  field.  Hence  the  nitrogen  spin  is  not  quantized  along  the 


5 


B  axis  and  any  calculated  rigid  lattice  second  moment  which  assumes  Bq 
quantization  will  be  too  low  (9).  Secondly,  molecular  motion  will  lead 
to  rapid  relaxation  and  this  relaxation  will  motionally  average  and 
so  diminish  the  N-H  interaction.  Finally,  the  N-H  x-ray  distances  are 
not  very  reliable.  With^these  caveats,  we  therefore  present  in  Table  I 
naive  estimates  of  the  n  contribution  to  the  proton  rigid  lattice  second 

moment:  the  indicated  errors  reflect  the  quoted  uncertainty  in  the  x-ray 

14  1  2 

distances  only.  The  average  N-  H  contribution  is  4.3+1G  ,  greater 

than  the  difference  between  the  observed  extra-pair  moment  (3.1G  )  and 

2 

the  calculated  1.50G  extra-pair  proton-proton  contribution.  This  result 
may  also  suggest  the  x-ray  N-H  distances  are  somewhat  short. 


More  information  can  be  extracted.  If  the  three  proton  pair 
separations  were  markedly  different  as  indicated  by  the  x-ray  results, 
then  the  resulting  proton  spectrum  would  correspond  to  three  Pake 
doublets  of  different  widths.  The  zero  crossing  technique  measures 
some  intermediate  pair  separation;  moreover,  when  the  extra-pair 
spectral  broadening  is  determined  by  dividing  out  the  intra-pair  fid 
of  Eq.  [2],  then  any  deviations  from  a  unique  Pake  splitting  will 
appear  as  an  apparent  extraneous  broadening.  We  can  roughly  estimate 
the  second  moment  of  apparent  broadening  by  the  standard  deviation  of 
the  three  intra-pair  moments  inferred  f^om  the  x-ray  distances.  The 
standard  deviation  is,  from  Table  I,  3G  ,  a  value  which  is  too  large 
to  account  for  the  difference  between  the  observed  extra  pair  moment 

2  2 
(3.1G  )  and  the  calculated  proton-proton  extra-pair  moment  (1.50G  ). 

Thus  these  NMR  second  moment  results  require  somewhat  larger 
and  more  nearly  equal  proton  pair  separations  than  those  determined 
by  x-ray.  Unfortunately  we  cannot  tell  with  these  data  whether  the 
difference  between  the  observed  extra-pair  second  moment  and  the 
calculated  proton-proton  extra-pair  moment  arises  primarily  from 

14  1 

different  Pake  splittings  or  from  the  N-  H  contribution  to  the 
second  moment. 


The  reconciliation  of  the  x-ray  and  NMR  derived  proton  separations 
is  not  so  central  to  this  report.  More  importantly,  as  we  have 
already  shown,  the  temperature  independence  of  both  the  extra-pair 
and  intra-pair  second  moments  rules  out  from  -60  to  230  C  any  molecular 
motion  which  might  modulate  the  dipolar  couplings  at  a  rate  greater  or 
equal  to  these  dipolar  field  frequencies.  Slower  motion  proceeds,  however, 
as  shown  in  Section  III  .  C. 

B.  Zeeman  Spin-lattice  Relaxation  (T^) 


The  Zeeman  spin-lattice  (longitudinal)  relaxation  time  is 
determined  by  molecular  motions  at  the  Larmor  frequency,  in  this  case 
60  MHz.  For  TATB,  proton  T^  measurements  were  performed  over  -80 

to  210°C.  In  addition,  a  separate  TATB  sample  was  pumped,  pressed  and 
then  y-irradiated  with  a  3  MRad  dosage;  T^  was  determined  for  this 
sample  at  high  and  low  temperature. 
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At  113°C  the  magnetization  data  are  shown  as  a  semilogarithmic 

plot  in  Fig.  3  and  as  magnetization  versus  t^^  in  Fig.  4.  Here  the 
magnetization  is  measured  by  sampling  the  height  of  the  initial  portion  of 
the  free  induction  decay:  this  is  proportional  to  the  integrated 
magnetization  across  the  entire  proton  spectrum.  The  decay  is  clearly 
nonexponential.  For  nuclear  relaxation  induced  by  rapidly  fluctuating 
electron-nuclear  dipolar  field  of  uniformly  distributed  paramagnetij 
centers,  one  naively  expects  the  magnetization  to  recover  as  (t/T^) 

for  short  times  (10) .  Of  course,  other  mechanisms  and  in  particular 
a  nonuniform  distribution  of  paramagnetic  centers  can  lead  to  non¬ 
exponential  recovery;  also,  over  a  sufficiently  limited  range  of 
magnetization  recovery,  almost  any  recovery  law  can  be  fit  (11) . 
Nevertheless,  in  Fig.  4  the  relation  M  oe  (t/T^)  ,  with  T^=35  s  ,  better 

represents  the  data  than  does  an  exponential  recovery. 

As  a  test  to  see  if  there  were  a  distribution  of  relaxation  times 
across  the  TATB  spectrum,  we  examined  further  these  113  C  data  as  well 
as  the  y-irradiated  sample  at  211  C.  The  fids  were  Fourier  transformed 
and  magnetization  plotted  versus  the  square  root  of  time  for  the 
center  of  the  spectrum,  one  peak,  and  one  shoulder;  in  the  Pake 
doublet  the  shoulder  is  twice  as  far  off  resonance  as  the  peak. 

The  center  of  the  pattern  corresponds  to  proton  pairs  aligned  at  54.7° 
relative  to  the  static  field  B  .  The  shoulder  corresponds  to 
orientations  parallel  to  the  field,  while  the  peak  is  a  mixture 
and  35.3  orientations.  The  decays  were  reasonably  fit  by  a  (t) 
dependence  even  out  to  the  point  that  the  fractional  magnetization  had 
recovered  by  about  one  half.  The  T^  values  corresponding  to  the  one 
half  point  (Table  II)  show  a  maximum  of  a  20%  variation  across  the  line: 
for  both  specimens,  the  fastest  relaxing  part  of  the  line  corresponds 
to  proton  pairs  parallel  to  the  field  Bq.  The  longitudinal  relaxation 

time  of  the  irradiated  sample  has  been  reduced  by  over  a  factor  of  two 
from  the  unirradiated  one.  It  is  more  convenient  to  characterize  the 
longitudinal  relaxation  data  at  each  temperature  by  a  single  parameter 
rather  than  record  the  slight  T^  dependence  across  the  spectrum  and  so 
we  report  the  T^  data  in  Figure  5  as  the  time  for  the  fid  magnetization 
to  recover  to  one  half  of  its  original  value. 

The  T^  data  are  essentially  temperature  Independent,  though  there 
is  perhaps  a  slight  trend  towards  longer  relaxation  times  at  higher  tempera¬ 
tures.  These  long  relaxation  times,  the  recovery  of  the  magnetization 
with  nearly  a  one  half  power  of  time  and  the  decrease  in  T^  on  y- 
irradiation  all  indicate  that  paramagnetic  centers  rapidly  relax 
nearby  protons  nuclear  spins;  by  mutual  spin  flips  the  proton  magneti¬ 
zation  diffuses  away  from  these  centers  to  build  up  macroscopic 
magnetization.  The  slight  T^  temperature  dependence  may  represent  less 
efficient  paramagnetic  relaxation  at  higher  temperature  or  a  slight 
reduction  in  the  spin  diffusion  rate  which  is  proportional  to  T^  ^ . 

The  virtually  temperature  independent  T^  values  show  that  there  is 

no  substantial  molecular  motion  at  high  frequencies  (60  MHz).  We  turn 
now  to  an  experiment  sensitive  to  slower  motions. 
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C.  Dipolar  Spin-lattice  Relaxation  (T^) 


The  rate  at  which  a  state  of  dipolar  order  relaxes  back  to  the 
lattice  temperature  can  provide  measure  of  slow  and  ultra  slow  mole¬ 
cular  motions  (12) .  The  decay  of  dipolar  magnetization  in  TATB  is 
generally  nonexponential  and  we  characterize  it  by  the  time  to  decay 
to  the  one  half  value;  however,  above  about  150°C,  the  decay  becomes 
more  and  more  exponential  and  for  those  temperatures  we  quote  also 
the  1/e  time.  As  shown  in  Fig.  5,  T.  slowly  diminishes  by  about  a 

o  o 

factor  of  four  from  -83  to  115  C  and  then  drops  rapidly  above  150  C. 

The  shortest  T^D  (1/e)  we  observed  is  3.7  ms  at  205°C. 

Dipolar  order  is  relaxed  by  any  motion  which  alters  the 
instantaneous  local  dipolar  field  seen  by  a  nuclear  spin.  For  vacancy 
diffusion,  the  general  relation  is  (12) 

tid  ■  Is 

2(l-pd)  [4] 

where  is  the  molecular  jump  time  and  (1-p^)  measures,  roughly 

speaking,  the  fractional  change  in  the  square  of  the  local  field  due 
to  that  diffusional  jump.  For  example,  p,  =  0.266  for  vacancy 

Q 

diffusion  to  nearest  neighbor  sites  in  a  powder  of  body-centered  cubic 
crystallites  (12). 

For  rotation  of  a  spin  pair,  the  relaxation  due  only  to  the 
modulation  of  the  pair  interaction  is 


T  =  T 

ID  _r  [5] 

d-Pr) 

For  three-fold  reorientation  perpendicular  to  the  internuclear  spin 
vector,  the  value  pr  =  -5/22  holds  (12).  This  analysis  would  apply 

for  TAT?  ^orientation  about  the  three-fold  axis. 

Pair  exchange  (two-fold  reorientation)  will  also  modulate  the 
local  field  arising  from  nuclei  outside  the  pair,  although  the 
efficiency  of  this  relaxation  process  is  reduced.  We  estimate  the 
efficiency  (1-p^)  of  this  process  to  be  equal  to  the  ratio  of  the 

proton-proton  extra  pair  contribution  to  the  second  moment  (1.5G^ 

2 

theoretical  value)  divided  by  the  total  proton  second  moment  (9  G 
experimental  intra  pair  moment  plus  the  1.5  G  ).  Hence  1-p^  21  1/7  for 
the-extra  pair  contribution  to  dipolar  spin-lattice  relaxation. 

Thus  the  difference  between  very  efficient  three-fold  axis 
reorientation  of  TATB  (1-p^  ”  27/22)  and  inefficient  180°  amino  proton 
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flips  (1-p^  -  1/7)  is  less  than  an  order  of  magnitude.  So,  indepen¬ 
dently  of  the  details  of  the  motion,  we  may  set  1-p  *»  1  in  order  to 
use  T^d  as  a  direct  measure  of  the  molecular  correlation  time  T  of 
some  unspecified  molecular  process: 


Tq  exp  H/kT 


[6] 


From  a  fit  of  the  T^D  (1/e  values)  at  the  four  highest  tempera¬ 
tures  on  Fig.  5,  we  find:  Tq  *  2x10  and  H  ■  120  kJ/mol.  The 
significance  of  this  activation  energy  is  discussed  in  Section  E. 

D.  Rotating  Frame  Spin-lattice  Relaxation  (T^  ) 


In  Fig. 

for  an  rf  field  strength  of  50  kHz. 
compressed  TATB  specimen  degraded  at  about 


5  we  also  show  the  limited  data  set  for  the  TATB  proton 

As  mentioned  earlier,  the 
l°C  (corresponding  T, 

data  are  not  shown)  and  damaged  the  probe,  so  precluding  further 
experiments. 


IP 


The  decay  of  rotating  frame  proton  magnetization  is  found  to  be 
nonexponential  in  TATB  and  the  T^  value  is  quoted  as  the  time 

for  the  magnetization  to  decay  to  one  half  of  its  initial  value.  The 
T^p  values  roughly  parallel  the  behavior  of  T^D  and  there  is  a  similar 
drop  at  high  temperature.  In  Fig.  5  the  T^D  temperature  dependence  has 
been  drawn  (dashed  curve)  through  the  few  T^p  data;  the  reasonable  fit 
shows  that  the  T^p  relaxation  process  is  also  consistent  with  the 
activation  enthalpy  of  120  kJ/mol. 


The  magnitude  of  the  T^p  values  is  not  fully  understood.  Experi¬ 
mentally,  T.  /T^Q  Z.  20  at  high  temperature  (T>170°C,  1000/T  <_  2.25) 
while  for  the  two  T^  data  at  lower  temperature,  T^/T^  *  2.5-3.  For 

illustration  only,  we  consider  isotropic  reorientation,  which  is  the 
most  efficient  dipolar  relaxation  process. 

In  a  powder  of  homonuclear  spins  irradiated  on  resonance  (13) , 


T 

C 


M„ 


2  2 
1+40).  x 
1  c 


[7] 


Here  is  the  second  moment  modulated  by  the  motion,  tc  the  correlation 

time  and  co^  is  the  rf  field  strength  in  rad/sec.  For  slow  motion, 

co.  T  »1  and  with  T,_  ~  T  ,  then 
1  c  ID  —  c 
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For  a  total  proton  second  iroaent  of  about  12  G  (experimental  values 

2  2 

of  9  G  intra-pair  and  3.1G  extra-pair),  Eq.  [8]  predicts  T^  /T^  *  45, 

a  factor  of  two  greater  than  that  observed.  Isotropic  reorientation 
is  unlikely  to  occur  in  TATB.  The  more  restricted  motions  are  less 
efficient  at  producing  T^  relaxation  and  also  generally  less  efficient 

at  producing  T^D  relaxation.  A  more  careful  calculation  considering  the 
detailed  dynamics  of  possible  relaxation  mechanisms  would  be  required  to 
extract  more  information  from  the  experimentally  determined  ratio  /T^. 

TATB  is  not  a  simple  spin  system,  involving  as  it  does  strongly  coup?ed^ 

possibly  inequivalent  proton  pairs  further  coupled  to  the  quadrupolar  N 
nuclei  as  well  as  to  the  rest  of  the  lattice.  For  now  we  note  only  that  at 
high  temperature  the  observed  ratio  is  broadly  consistent  with  the 
determination  of  T^  and  T^D  by  a  single  relaxation  time,  described  by  Eq.  [6]. 

At  lower  temperatures  below,  say,  150°C,  T^  is  only  a  few  times 
larger  than  T^  .  Unfortunately  from  those  data  we  cannot  yet  determine 
the  T^  or  T1D  relaxation  mechanism  in  this  region. 

E.  Molecular  Motion  in  TATB 

From  the  measurements  of  T^Q  at  high  temperature,  where  molecular 

motion  dominates  the  relaxation,  we  have  estimated  the  jump  time  for 
the  motional  process  to  be 

T  =  2xl0-16  exp  (120  000/RT)  [9] 

where  R  is  in  units  of  J  mol  ^  deg  ^  and  x  is  in  seconds.  The  TATB 

molecule  has  various  possibilities  for  molecular  motion:  translational 

diffusion,  molecular  rotation  about  the  three- fold  symmetry  axis,  or 

even  rotation  of  nitro  as  well  as  amino  groups.  For  diffusion,  the 

activation  enthalpy  for  planar  molecules  is  generally  somewhat  greater 

than  twice  the  latent  heat  of  sublimation,  as  exemplified  by  the 

entries  in  Table  III  for  benzene  and  naphthalene  (14) .  A  reliable 

value  for  the  latent  heat  of  sublimation,  AH  ,  ,  of  TATB  does  not  exist 

sub 

in  the  literature  (15) .  Thus  to  estimate  the  activation  enthalpy  for 

translational  diffusion  in  TATB  AH  ,  itself  must  be  estimated,  as  is  done 

sub 

below  in  three  ways.  First,  TATB  is  to  a  good  approximation  the  same 

size  as  coronene,  for  which  experimental  and  theoretical  estimates  of 

AH  ,  exist  (Table  III),  about  140  kJ/mol.  (16).  Second,  from  mass 
su  b 

spectroscopic  measurements,  it  has  been  found  (15)  that  TATB  begins  to 

show  an  appreciable  vapor  pressure  (about  10  ^  torr)  at  200°C;  use  of 
this  datum  with  the  approximate  equation  (17) 


[10] 


log1(J  p(torr)  *  13  -  AHgub/(2. 303RT) 

yields  a  value  of  145  kJ/mol  for  AHgub>  Finally,  a  table  (18)  of 

functional  group  increments  for  AH  ,  is  used  to  estimate  AH  ,  of 

sub  sub 

180  kJ/mol.  Thus  the  activation  enthalpy  for  diffusion,  expected  to  be 

about  twice  AHgub,  should  be  about  300  kJ/mol .  Since  this  is  much 

larger  than  theactivation  enthalpy  of  120  kJ/mol  found  in  this  study, 

the  translational  diffusion  process  is  eliminated  as  the  source  of  the 

high  temperature  relaxation  effect. 

In  Table  III,  activation  enthalpies  for  molecular  rotation,  AH  , 

rot 

are  also  gathered.  Coronene,  a  molecule  with  a  six-fold  axis  of 
symmetry,  requires  only  25  kJ/mol  for  rotation  (16),  while  naphthalene, 
with  only  a  two-fold  axis,  requires  105  kJ/mol.  (16) .  For  neither  of 
these  do  hydrogen  bonding  or  molecular  electric  dipole  moments  play  any 
role  in  the  potential  energy  function.  TATB,  with  a  three-fold  axis, 
and  roughly  the  same  moment  of  inertia  as  coronene,  should  resemble 
coronene,  except  for  the  possible  contributions  of  polar  groups  and 
hydrogen  bonding.  Thus  25  kJ/mol  is  a  reasonable  lower  bound  to  the 
enthalpy  of  activation  for  rotation  in  TATB.  But  a  rotation  of  the  amino 
group  by  180°  is  also  a  possible  contributor  to  relaxation,  though  much 
less  effective. 

A  strategy  to  discriminate  between  the  two  rotational  possibilities 

is  a  study  of  relaxation  in  the  rotating  reference  frame  as  a  function 

of  rf  field  strength,  which  we  have  begun  but  not  yet  completed. 

Further  studies  at  higher  temperature  may  be  useful,  although  limited. 

Though  TATB  does  not  melt  until  about  440°C  (19) ,  the  degradation  at 

350  C  (20) ,  as  well  as  the  degradation  we  observed  at  around  300°C  in 

the  compressed  sample,  provide  a  very  inconvenient  ceiling  on  the 

accessible  temperatures.  By  extrapolation  of  our  data,  we  should 

expect  at  300  C  that  substantial  line  narrowing  will  occur  as  the 

molecular  jump  time  becomes  comparable  to  the  T2  (zero  crossing)  time. 

Further  temperature  dependent  rotating  frame  relaxation  time  results 

would  be  helpful  in  assessing  the  details  of  the  motion,  although 

we  estimate  that  the  T,  minimum  would  not  occur  until  about  365  C  for 

lp 

the  rf  field  strength  of  50  kHz.  (The  observation  of  the  T^  minimum 
is  strictly  precluded,  as  it  probably  occurs  at  around  575°C.)  At 
lower  temperature,  nuclear  spin  echo  as  well  as  esr  studies  could  more 
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fully  identify  the  role  of  N  relaxation  and  paramagnetic  impurity 
relaxation. 

IV.  CONCLUSIONS 

Below  about  150°C  there  is  no  molecular  motion  in  TATB  which  is 
detectable  by  these  NMR  experiments.  The  determination  of  the  dipolar 
spin- lattice  relaxation  times  (T.^),  and  rotating  frame  relaxation 


times  (T^p) ,  which  are  extremely  sensitive  to  slow  motions,  has  pointed 

to  a  motion  with  an  activation  enthalpy  of  120  kJ/mol.  The  correlation 
time  x  for  this  motion  can  be  extracted  from  the  T1D  data  only  on  the 

assumption  of  a  particular  motional  model;  nonetheless,  it  is  a  fair 
estimate  to  set  x  -  T,_.  and  from  this  we  inferred  a  prefactor  in  the 

-16 

Arrhenius  relation  of  2x10  s.  At  room  temperature  this  motion, 
responsible  for  the  relaxation  at  high  temperature,  corresponds  to 

jumps  occurring  approximately  every  five  days. 

The  motion  responsible  for  the  high  temperature  (above  150°C)  T^Q 
and  T^p  relaxation  cannot  be  unambiguously  determined;  180°  flips  of 

amino  protons  and  three-fold  axis  reorientation  are  possibilities.  On 
comparison  with  other  characterized  organic  molecules,  the  activation 
energy  we  observe  is  smaller,  by  more  than  a  factor  of  two,  than  that 
expected  for  translational  diffusion  and  so  the  relaxation  data  we 
present  is  inconsistent  with  translational  diffusion  at  these 
temperatures . 

The  T^  values  and  Pake  doublet  lineshape  are  characteristic  of 

strongly  coupled  dipolar  pairs  which  are  only  weakly  linked  to  the 
other  nuclei  in  the  crystal.  As  an  indication  of  hydrogen  bonding, 
the  proton  spectrum  is  shifted  downfield  to  6  »  34  +  8  ppm.  On  the 
assumption  that  the  three  amino  proton  pair  separations  are  identical, 
we  infer  from  the  T2  data  a  proton  separation  of  0.185  +  0.001  nm,  a  value 
outside  the  x-ray  results  of  0.1553,  0.1577  and  0.1726  +  0.02  nm;  we  feel 
these  latter  separations  are  systematically  in  error. 

The  broadening  of  the  Pake  pattern  can  be  described  by  a  Gaussian 
2  2 

with  a  3  G  second  moment.  Of  this,  1.5  G  arises  from  proton-proton 
coupling  from  outside  the  amino  pair.  Nitrogen-proton  dipolar  coupling 

2 

may  account  for  4.3  G  and  a  distribution  of  amino  proton  pair  separations, 
as  given  by  the  x-ray  results,  could  also  account  for  3  G  .  The  experi¬ 
mentally  observed  broadening  indicates  that  the  amino  proton  pair 
separations  are  not  so  disparate.  Both  the  intra  and  extrg-pair  proton 
second  moments  are  temperature  independent  over  -80  to  230  C. 

The  Zeeman  spin-lattice  relaxation  times  T^  are  also  essentially 
temperature  independent  in  the  range  accessible  and  are  most  likely 
determined  by  the  relaxation  of  paramagnetic  centers  within  the  crystal. 
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f  (kHz)  — ► 

Fig.  1  —  Proton  NMR  spectrum  of  TATB  (a)  as  received  and  (b)  after  pumping  at  200°C.  The  sharp  peak  in 
(a)  is  ascribed  to  water  present  in  the  processed  material.  A  rather  narrow  proton  resonance  from  the  spectro¬ 
meter  probe  has  been  subtracted  off  both  (a)  and  (b),  but  the  cancellation  is  not  perfect  and  the  features  in  the 
valley  at  about  75  kHz  offset  arise  from  this  probe  resonance  and  not  TATB.  The  Pake  doublet  is  characteris¬ 
tic  of  strongly  coupled  isolated  spin  pain,  here  of  the  amino  protons.  The  splitting  of  28.5  kHz  is  derived  from 
experimentally  determined  time  of  the  first  zero  crossing  of  free  induction  decay  (fid)  and  represents  the  sepa¬ 
ration  of  the  sharp  peaks  in  the  Pake  doublet. 


second  moment.  Both  second  moments  are  essentially  temperature  independent. 


Fig.  3  —  Magnetization  versus  time  data  for  Zeeman  spin-lattice  relaxation  (Tj)  at  113°C.  The  recovery 
of  magnetization  is  nonexponential;  see  also  Figure  4.  The  total  magnetization  is  represented  by  using 
the  initial  height  of  the  fid. 


Vt  Usee)  ■  ■-* 

Fig,  4  —  The  Tj  magnetization  data  of  Figure  3  are  better  fit  by  an  M  versus  (t/Tj)1/j  relation.  The 
solid  line  corresponds  to  Tj  -  35  s. 
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Fig.  5  —  This  figure  summarizes  the  four  relaxation  times  measured:  T2  (by  zero  crossing),  dipolar  spin-lattice  time 
Tjq,  rotating  frame  spin  lattice  time  Tjp  and  Zeeman  spin-lattice  time  Tj.  The  nonexponential  decays  for  T^q, 

T  jp  and  Tj  are  characterized  by  the  decay  time  to  the  one  half  point,  except  at  high  temperature  for  T^q  where  the 
(1/e)  values  are  also  shown.  The  only  evidence  of  molecular  motion  comes  from  the  sharp  drop  in  T^q  and  T^„ 
above  about  150°C:  the  activation  enthalpy  is  120  kJ/mol.  A  7-irradiated  specimen  shows  a  reduced  by  a  factor 
of  two,  so  indicating  paramagnetic  relaxation  as  the  Tj  mechanism. 
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the  unit  cell:  here  they 


TABLE  II 


Relaxation 

Across  the  TATB 

Proton  Spectrum 

Temperature 

113°C 

211°C 

(Y-irradiated) 

T^Cs) 

T^Cs) 

f  idb 

36 

16 

on  resonance 

32.5 

15 

peak 

36 

15 

shoulder 

41.5 

18 

is  shown  here  as  the  time  to  recover  to  one  half  of 
the  value  of  the  magnetization  at  thermal  equilibrium, 
corresponding  to  a  pulse  delay  time  of  1000s. 

bThe  initial  part  of  the  free  induction  decay  (fid) 
measures  the  integrated  spectrum. 
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TABLE  III 


Si 

Enthalpies  of  Activation  and  Sublimation 
(in  kJ/mol) 


Process 

Rotation 

Sublimation 

Diffusion 

Substance 

AH  „ 
rot 

(expt)  AH  (calc) 

rot 

AH  k<expt) 
sub 

AH  .<CalC> 
sub 

AHdiff. 

Benzene 

16 

8.8-13 

44.8 

38-44 

96b 

Naphthalene 

105 

65-88 

67.4-82.0 

63-71 

179C 

Coronene 

25 

8-12 

144 

138-151 

- 

TATB 

_ 

_ 

145d 

180e 

— 

aData  taken  from  Ref.  16  unless  indicated  by  other  superscript. 

bR.  Fox  and  J.  N.  Sherwood,  Trans.  Farad  Soc.  67 ,  3364  (1971)  and 
F.  Noack,  M.  Weithase  and  J.  von  Schutz,  z.  Naturf.  30a,  1707  (1975). 

CJ.  N.  Sherwood  and  D.  J.  White,  Phil.  Mag.  1£,  975  (1967). 

^Deduced  from  vapor  pressure  estimate;  see  text. 

Calculated  from  functional  group  increments  as  described  in  text. 
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